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Abstract 

The primary observable in dark matter direct detection is the spectrum of scattering events. We 
simulate multiple positive direct detection signals (on germanium, xenon, and argon targets) 
to explore the extent to which the underlying particle physics, manifested in the momentum 
dependence of the operator mediating the scattering, can be extracted. Taking into account 
realization (Poisson) noise, a single target nucleus with 300 events has limited power to discriminate 
operators with momentum dependence differing by q^^ for a wide range of dark matter masses from 
10 GeV to 1 TeV. With the inclusion of multiple targets (or a factor of several more events on a 
single target), the discrimination of operators with different momentum dependence becomes very 
strong at the 95% C.L. for dark matter candidates of mass 50 GeV and above. On the other hand, 
operator discrimination remains poor for 10 GeV candidates until multiple experiments each collect 
1000 or more events. 
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I. INTRODUCTION 

Dark matter (DM) dominates the dynamics of matter on large scales in our universe, 
yet the theory which describes its underlying interactions, both with itself and with the 
Standard Model (SM) sector, remains hidden. Viable, testable models that are capable of 
both describing the DM and satisfying the observed constraints abound in the literature. 
The most popular DM candidate is the weakly interacting massive particle (WIMP). A 
well-studied example of a WIMP is the lightest supersymmetric particle in the minimal 
supersymmetric standard model (see e.g. PP for a review). This candidate satisfies the 
observational constraints that it is long-lived, cold, and weakly interacting. Its mass is 
typically predicted to be between 10 GeV and 1 TeV, and its interaction with ordinary 
particles is mediated by an operator that gives rise to momentum-independent interactions. 

Direct detection experiments, such as XENONIOO [2] and CDMS-II |3j, have achieved a 
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high level of sensitivity to a WIMP in this mass range. In general, the focus has been on 
ruling out DM candidates with the most standard characteristics - momentum independent 
scattering in the mass range well above 10 GeV. However, some recent experimental hints 
have led to more detailed investigations of other possibilities. In particular, the DAMA- 
LIBRA [1], CoGeNT [5], and CRESST-II [6] experiments have excesses which may be 
consistent with a DM mass of order 10 GeV, leading to a renewed theoretical interest in 
DM candidates with mass in this range. However, these results are in tension with the null 
results of XENONIO [3 [8], XENONIOO [2], and CDMS [9]. This has revitalized interest in 
DM candidates that have non-standard features that might allow one to evade the constraints 
from the null experiments. In particular, operators that mediate interactions that give rise 
to momentum-dependent scattering rates, first considered in [T0HT2] , have received renewed 
interest [IMH]. 

Independent of these results, direct detection experiments with increasing sensitivity are 
continuing to probe the WIMP DM hypothesis. The XENON collaboration will soon begin 
construction on a XENONIT phase [191; ^ 3.6 ton liquid argon detector operated by the 
DEAP/CLEAN collaboration, DEAP-3600 [20|, is under construction; and 100 kg and 1 
ton cryogenic germanium detectors are being planned [21] by the Super-CDMS/GEODM 
collaborations. If there is a positive signal from one of these experiments, the signal must 
be confirmed by multiple experiments with different targets. Furthermore, one would like 
to determine experimentally whether the scattering is in fact mediated by a standard spin- 
independent operator or whether the underlying particle physics is more complex, with non- 
standard types of momentum dependence in the scattering. Additionally, recent theoretical 
work has focused on the feasibility of extracting different types of information from scattering 
events [221123]. 

In this paper we explore the dark matter inverse problem; that is, the capability of direct 
detection experiments to extract the underlying particle physics mediating the scattering of 
a DM particle. In practice this question boils down to how well the momentum dependence 
in an operator can be mimicked by other operators for a given target and DM mass. We 
simulate detection events using different types of interactions with varying momentum and 
velocity dependence, then we fit these events by a wide variety of interactions. In principle, 
one can distinguish interaction types by just checking the recoil spectra for an individual 
target. For example, if DM scatters with the target nucleus through an interaction with 
a positive momentum dependence, one would expect events to drop in low recoil energy 
bins, while the event rate increases exponentially towards low recoil energy for a standard 
momentum independent interaction. However, if only one target is available, the capability 
for distinguishing operators through the spectra in this way is limited by the experimental 
noise and detector threshold. With multiple nuclear targets one may also compare the overlap 
of the DM preferred regions on different target nuclei to determine the correct operator. We 
find that, for similar cross-sections, analysis of high mass DM particles leaves distinctive 
qualitative imprints that allow one to extract the momentum dependence in the scattering. 
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When two operators cannot be distinguished between each other, it is because they have a 
very similar momentum dependence. On the other hand, operator discrimination for a 10 
GeV candidate is poor. We conclude that a significantly lower threshold, or much improved 
statistics, than is available for the next generation of experiments will be necessary in order 
to extract the particle physics mediating the scattering for a light DM candidate. 

We begin the paper by briefiy reviewing direct detection basics. In Section III we discuss 
our methodology, including detector effects, analysis methods, and conventions. In Section 
IV we go over our results. Finally, we summarize our conclusions in Section V. 



II. PRELIMINARIES 



We compare the spectra for standard spin- and momentum-independent scattering against 
operators with momentum dependence. We call operators with n additional powers of 
momentum transfer g^-type scattering, where n = ±2, ib4. We also consider scattering 
via anapole and dipole operators, which feature mixed momentum dependence (constant 
and dependence for the former operator, and and dependence for the latter). In 
this section we review the rate relations for the various types of operators and define our 
conventions for reporting results. 

For generic scattering the observed differential rate of observation of DM particles may 
be written 

dR / Po f r/- ^ \ /7-iN d(^ DM 



dE 



R 



(Er) = — [ £v V fivo,Ve;V,nin,Vesc)eesiER)^^^^^^^{v, A, Er). (1) 

mDurriN J dER 



Here, po is the local DM density; moM is the DM mass; niN is the detector nucleus mass; the 
velocity distribution is given by /(fo, Ve] Vmin, Vesc)', ^es{ER) is the detector efficiency, which 
may depend on the energy of recoil; and the differential cross-section daDu-N / dER describes 
the interaction. The total number of events is given by integrating the differential rate as 

T-imax.i 

/■ « dR 

iV,(E— )=/ dER—e,^^, (2) 

Jri 

where eexp is the exposure, generally given in terms of kilograms of exposed material multi- 
plied by the days of duration of the experiment. 

All of the particle physics information present in Eq. ([T]) is contained in the final term, 
da^yi-N / dER, while the astrophysical considerations are refiected in po and /(fo, Ve, Vmin, Vesc)- 
The remaining effects, such as ees{ER), eexp, and m^, are detector-specific. The objective of 
this paper is to extract information about the operator content of da^yi-N / dER based solely 
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on observables. We fix the astrophysical parameters by the best current observations and 
assume the detector effects are well understood. In particular, we take po = 0.4 GeV/cm^ 
and assume a Maxwell-Boltzmann velocity distribution with a smooth exponential cutoff at 
the escape velocity: 



esc 



f{VQ,Ve]Vmin,Vesc)={ ^ ^ , (3) 

) V > Vesc 



where we take v^sc = 544 km/s, Vq = 220 km/s, and = 232 km/s. The normalization for 
this distribution can be found in the appendix of [TB] . We find that our results are insensitive 
to changes in these velocity parameters. Note here that we are only interested in extracting 
the correct particle physics interaction. This differs from recent work where the goal was to 
marginalize over astrophysical uncertainties while extracting the basic physical parameters 
[25] . Astrophysical uncertainty will enlarge the preferred region in the relevant parameter 
space for a single experiment, but our results involve a comparison of multiple experiments 
for which the astrophysics are the same. We have checked that our results are unchanged 
as long as we use the same inputs for the different experiments and do not marginalize over 
the unknown quantities. 

The differential scattering cross-section is related to the DM-nucleus scattering cross- 
section (Tat via 

where /xat is the reduced mass of the DM-nucleus pair and v is their relative velocity. For 
standard spin- and momentum-independent scattering, which we denote "std", the results 
are quoted in terms of the cross-section for scattering off a nucleon, a^'^: 

(^N =^n TY 7-2 FAAEr), (5) 

where /i„ is the reduced mass of the DM-nucleon pair. We take fp = fn = 1 here, and 
assume a standard Helm form factor, Fi{A,Eji) = 3ji(gro)/(5'ro)e~^^^^^^™ with s = 0.9 
and ro = ((1.23^1/3 - 0.6)^ + 7/3(0.527r)2 - 5s^Y^^ fm 

For spin- independent, momentum-dependent scattering, the results will be reported ac- 
cording to the convention |13j : 

(!)"<. (6) 

In this work, we take go = 50 MeV for all values of We allow (T„ to float for both 



^ Particle physics realizations of momentum-dependent scattering cross sections are discussed in jl3j. Note 
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the standard and type operators to satisfy normalization conditions described below. In 
general, this will lead to (T„ ~ (9(10"'^^) cm^. Note that the momentum-transfer dependence 
of the g^-type operators serves to introduce extra powers of the recoil energy, not extra 
powers of the DM velocity. 

Lastly, we define our conventions for operators with mixed momentum dependence, the 
anapole (xT^TsX^m) and dipole (x'^'^'^X^fiu) operators. The scattering cross-sections for 
these operators are [16] 



dip _ 



4q\^Z^F^{A,ER 



(7) 



Here, M is the mass of a mediator particle that couples an off-shell photon to the DM 
particle, while A is a confinement or compositeness scale that describes the magnetic dipole 
moment physics; we take it to be 100 MeV in all of our analysis. The mediator mass M 
is a priori unknown. In general, it floats around a fiducial value of a few hundred GeV. 
J is the spin of the target nucleus, is the magnetic moment of the target nucleus, and 
bn = e/2mp is the Bohr magneton. The form factor F2{A,Eii) for the spin coupling term 
is important for spin-odd nuclei and is taken from [26]. Since the isoscalar and isovector 
couplings are model-dependent, we choose to take Op = a„ in this work. We find that the 
results are qualitatively insensitive to these exact values. 

The normalized interaction strength that we display on the plots for both anapole and 
dipole moment operators is 

^„ = ^ 1.05 X 10-39 (—X( . ' cm'. (8) 

47rM4 \mnj V 400 GeV / ^' 

The fiducial value of M and thus of a„ differs depending on the type of interaction and 
the mass of the DM, but these values are chosen to obey the same normalization conditions 
alluded to above. Our normalization conventions are described in greater detail in the next 
section. 



that |13) assumes that the mediator mass is heavier than the momentum transfer so that n = 2,4 are 
positive. If the mediator mass is smaUer than the momentum transfer, we can have n = —2, —4. 
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III. METHODOLOGY 



We simulate mock experimental ( "input" ) data, and compare the goodness-of-fit of theo- 
retical ("fit") spectra against the input data. Ideal experimental data is first generated via 
a model, then smeared according to both Poisson statistics and the finite energy resolution 
of the experiment. We then fit theoretical models (convoluted to take into account detector 
resolution) to the input data using a log-likelihood ratio test. We describe our methodology 
in this section. 



A. Detector and Statistical Effects 



In order to fully simulate a direct detection experiment, one must fold detector effects into 
the rate given in Eq. ([T]). A convolution integral (accounting for the finite energy resolution 
of the detector) and Poisson noise (expected in a counting type experiment) are the most 
important of these effects. For low mass DM in particular, these considerations are equally 
important in drawing conclusions. We describe here how we take these effects into account. 

Convolution. In order to compare a theoretical model against an input data set, the 
rates must be smeared to take into account detector effects and noise. Both the input 
data and the theoretical spectra are smeared according to the energy resolution, because a 
realistic detector does not have perfect energy resolution. The observed rate is well modeled 
by convolving the expected rate with an approximately Gaussian distribution whose width 
is a function of detector element and true recoil energy: 



dR 



dE 



R 



1 



dE 



R 



/271 adet{A,E'] 



exp 



{En - E'f ^ 



(9) 



For all isotopes present in a given detector we define the detector energy resolutions to be 



arfet(Ge, E) = ^(0.3)2 + (0.06)2E/keV keV, 



adei(Xe, E) = 0.6^/E/keV keV, 



adetiAr, E) = CTV^/keV keV. 



(10) 

(11) 
(12) 



In all cases we take a flat efficiency. The isotopic abundances are crucial for extracting 
information in the anapole and dipole cases. We list the spins, magnetic moments, and 
abundances of the important spin-odd nuclei for each element in Table [T], and take the 
abundances for the spin-even nuclei to be set by naturally occurring levels. 
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Element 


A 


Z 


J 


hN/hn 


Abundance 


Germanium 


73 


32 


9/2+ 


-0.879 


7.73% 


Xenon 


129 


54 


1/2+ 


-0.778 


26.4% 


Xenon 


131 


54 


3/2+ 


+0.692 


21.2% 



Table I. Physical properties of stable spin-odd nuclei present in detectors. The dominant argon 
nucleus has A = 40, and there are no spin-odd isotopes with appreciable abundance. 

Integrating the convoluted rate in Eq. (|9| over the recoil energy range in each bin i gives 
rise to binned observation numbers 



N^{E_ 



j-imax\ 
R ^^R ] 



dEi 



dR 



dE 



R 



-exp- 



(13) 



These binned observation numbers produce the observable quantities that make up both our 
input data and the data used in our fits. We will refer to the vector of data points generated 



by Eq. (13) as the "convoluted spectrum." 



Poisson Noise. We expect that our input data are one of an ensemble of many other 
data that could be produced by the expected model. Since the input data that come from 
any given experiment are produced by counting, we expect that the input data will deviate 



from the convoluted spectra given in Eq. ( 13 ) at a level determined by Poisson statistics. To 



simulate the observed events, we thus introduce stochastic Poisson noise distributed about 



the spectrum given in Eq. (13). 



We have checked that, as expected, the average CLC generated by many noisy data sets 



converges on a CLC described by data generated by inputs from Eq. (13) with no noise. 



Log-Likelihood Ratio Test. We use a log-likelihood ratio to provide an estimate of the fit 
parameters that best reproduce a set of input data. For input data rii in the i-th energy bin 
(which has been smeared in energy as described above and in the number of events according 
to Poisson statistics) and theoretical parameters 9 that produce a convoluted spectrum with 
values Ui, the log-likelihood ratio is 



TV 



i=l 



UiiO) — + rii In 



+ const.. 



(14) 



where if nj = the last term in the brackets is set to zeroj^ The parameters of interest are the 



DM-nucleon cross section and DM mass. The parameters that minimize Eq. ( 14 ) provide the 



^ The value of Vi should never be exactly zero because even at high energy the convolution integral will 
receive some contribution from events in the low energy tail of the convolution. 
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point of best fit, and contours of constant L values give confidence level regions. To define 
the 95% confidence level regions we take the region inside of which L < + ^-^95*^ °'^'; 
where for instance AL^5^°-^- = 5.99 pi] . 



B. Defining the Spectra 

We generate a vector of input data in the i energy bins, n^, for all three target nuclei 
assuming interactions derived from standard, anapole, dipole, g^, and operators. Recall 
that the input data are smeared in energy as well as with Poisson noise. These spectra 
are obtained for DM masses of 10, 50, and 250 GeV. This is then compared to the vector 
of theoretical data in the i energy bins, z/j, which are the convoluted spectra for standard, 
anapole, dipole, g^^, and g^"^ operators. 

For xenon and germanium experiments, we take the energy range of the experiment to 
be 5-60 keV, binned in a 1 keV interval from 5 to 13 keV, 2 keV from 13 to 25 keV, and 5 
keV from 25 to 60 keV. For the argon experiment, a higher threshold is assumed, taking the 
energy range of the experiment to be 20-60 keV, binned in 1 keV from 20 to 28 keV, 2 keV 
from 28 to 40 keV, and 5 keV from 40 to 60 keV. These bins are moderately coarser than 
current and upcoming experimental capabilities, but we find that the results are insensitive 
to bin size. We have checked that more than doubling the number of allowed bins and using 
bins as narrow as 0.4 keV does not affect our results. In all cases, properly accounting for 
the effect of the noise and the energy resolution are more important than binning choices for 
reaching robust conclusions. 

Our conventions for the total rate may be fixed in two ways. First, we take the exposure 
of all targets to be 2 ton-years and fix the total number of events on germanium to be 300 
for 50 and 250 GeV DM candidates and 100 for a 10 GeV DM candidate. We allow (T„ and 
an to vary so as to achieve the desired number of events on germanium. For a fixed exposure 
more events are observed on a xenon target due to the higher atomic number, while fewer 
are observed on an argon target. To check the robustness of our results and their sensitivity 
to statistics, we also calculate results when an equal number of events are obtained on all 
targets (300 events for 50 and 250 GeV DM test masses, and 100 events for 10 GeV DM test 
mass). This will allow us to disentangle the effect of statistics from the effects of particle 
physics. We find that while the constraining power of individual experiments changes, we 
are able to obtain good discrimination in both cases. 

Due to threshold effects, fewer events are expected to result for 10 GeV DM candidates 
for a given scattering cross-section. As a result of these fewer events, noise is an important 
limiting factor in our ability to definitively extract the correct operator mediating the 
interaction. For light DM, kinematics is also very limiting. For these reasons, our conclusions 
for light DM are less sensitive to the total number of events: even with as many as 1000 
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events on a single target (which for a germanium target corresponds to 20 ton-years of 
exposure for a cross-section (Jn ~ 10"*^^ cm^) we do not obtain good discrimination for a 10 
GeV candidate. 



IV. RESULTS 



There are two possible measures for operator discrimination that will be obvious from 
our results: the absolute goodness-of-fit per degree of freedom, Lmm/d.o.f., and the overlap 
of the 95% preferred regions from fits to spectra generated by different target nuclei. A 
mismatch of the input data and theoretical model should give high log-likelihood values, and 
comparing the Lmin/d.o.f. values for different trial operators can provide a clue to nature 
of the interaction that generated the data. Similarly, different operators will give rise to 
different preferred regions as observed on multiple targets. Overlapping confidence level 
curves (CLCs) can indicate that the momentum dependence scales with target nucleus as 
expected for the hypothesized interaction and disjoint CLCs can indicate the converse. As 
detailed below, considering the value of Lmin for the combined data from all targets gives a 
test that is sensitive to both of these measures. 

Two subsets of the results are shown in Figs. ([T]) and ^ for input data corresponding 
to standard and dipole mediated scattering. More complete results are available in the 
appendix and online [27]. The subsets shown in Figs. ([T]) and ^ are chosen to represent 
illustrative contrasts to the input data; we describe our particular choices in more detail 
below. Each figure corresponds to one set of input data (labeled by rii) generated according 
to the methodology described in the previous section. The individual panels of each figure 
correspond to fits of a theoretical model (labeled by ui) to the input data. We divide the 
plots in Figs. ([T]) and ^ so that the left panels show our results for the 10 and 50 GeV 
DM candidates for scattering on all three targets: xenon, germanium and argonj^ The right 
panels show the same information for 250 GeV candidates. In the appendix, all CLCs are 
shown on the same plot. The colors on the curves correspond to the minimum log-likelihood 
per degree of freedom|^for the fit of each operator to the data, with darker colors representing 
the better fits. We have checked that different realizations of noise do not change our ability 
to extract information about the particle physics. 

In this section, we will start by going through a simplified analysis of the scattering 
kinematics that explains the shapes of the CLCs in Figs. ([T]) and (|2]). Next, we will explain 
how the tests discussed here allow one to extract the correct particle physics that underlies 
direct detection events. We will do this both when all exposures are the same and when the 

^ We omit the argon curves in the 10 GeV case because argon events generated for the 10 GeV analysis were 
essentiaUy compatible with zero. Our argon curves are therefore closer to exclusion curves than preferred 
regions. In all cases the germanium and xenon CLCs fit safely underneath the argon "exclusion curve." 

^ For high mass, d.o.f.(Ge, Xe) = 19, while d.o.f.(Ar) = 16, corresponding to the number of bins minus 
the two fit parameters, ttidm and (T„. For low mass candidates we have many fewer degrees of freedom 
because fewer bins are filled: d.o.f.(Ge) = 6 and d]cQf.(Xe) = 1 — 5, depending on the operator. 



number of observed events is the same on all targets. This allows us to determine the effects 
of statistics on our results. Finally, we will compare models using the p- value test. This test 
will be a quantitative measure of both of the qualitative tests described above. 



A. Contour Shapes from Scattering Kinematics 



We discuss the shapes of the CLCs that we observe in Figs. ([T]) and (|2]). The kinematics 
of elastic scattering between DM particles and detector target nucleus provides a straight- 
forward handle on the shape of these preferred regions. For the purpose of illustration, we 
take standard scattering and simplify the differential event rate given in Eq. ([T]) with the 
limit Vf-sc oo, ignoring solar motion, terrestrial motion, detector efficiency, and detector 
resolution. This recovers the result pQ 

-{Er) ~ 2 — exp , , . (15) 



The balance between the exponential term and the linear term on the right-hand side of 



Eq. (15) gives most of the important kinematic information. This balance is the key to 
understanding the basic content of our plots. To better understand the implications of 
scattering kinematics on the operator discrimination, we consider two extreme limits first. 

If the DM mass is much less the target nucleus mass then fj,^ ~ ""^dm and we can 
approximate the differential rate by 

dR PqGn ( EnniN \ 



(IEr aA^^dm^o V '^'^dm'^o/ 

We can immediately see that this rate is very sensitive to the DM mass. If we increase the 
DM mass the scattering rate will increase sharply because of the exponential term. Our 
numerical results clearly refiect this sensitive dependence. From fits of the 10 GeV case 
shown in Figs. ([T]) and ^ we can see that all the 10 GeV contours are very narrow. By 
contrast, the event rate is linearly proportional to the scattering cross section and the shift 
has a rather mild dependence on the target mass. Hence, detectors have poor cross-section 
discrimination capabilities for the low mass case. 

At the other end of the mass range, if the DM is much heavier than the detector nucleus 
mass then the differential event rate is given to first order by 

dR ,„ ^ Po(Tn { Er \ 
; exn ^ 



dER y/Trm-DMm%Vo \ 2m ^v^ J 



(Er) - % exp . (17) 
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The exponential term is independent of the DM mass in this hmit, and it becomes an 
overall factor for any given large DM mass. The differential event rate now just scales as 
dR/dEji oc 0"„/mDM- In this limit the target nucleus is not able to determine the DM 
mass because one always can rescale (T„ to compensate the change of DM mass while keeping 
o'n /'mDM unchanged. This so-called high-mass degeneracy [23] is observed in Figs. ([T]) and 
The allowed mass range can be unbounded for all three targets when high-mass degeneracy 
is important. 

These features explain certain coarse aspects of the CLC plots. Now we examine the 
physics that leads to the separation of the CLCs in the parameter space, and we investigate 
what this can tell us about the particle physics underlying the data. We do this analysis 
first with equal exposures on all targets. We will also present results with the same number 
of events recorded by all targets. 



B. Operator Discrimination with Equal Exposures on All Targets 

For operator discrimination with equal exposures, the most important feature of Figs. ([T]) 
and ^ is that, when the input data and theoretical model match (nj = z/j), the CLCs overlap 
regardless of candidate mass. When the theoretical model is mismatched with the input data, 
the 50 GeV and 250 GeV CLCs separate, though much less separation is apparent in the 
10 GeV case. Thus for low mass candidates the overlap test does not strongly discriminate 
between models, while for high mass candidates the overlap test allows strong discrimination. 
Much of this power is derived from the strong statistics gained from scattering off of a xenon 
target. By contrast, we see in Fig. ([2]), where the momentum dependence of the operators 
differs by g^, that, for scattering off a single target nucleus, the log-likelihood per degree- 
of-freedom test can be limited by noise and does not provide a very good discriminant 
between theoretical models - the fit can be just as good on a single target with the incorrect 
theoretical model as it is with the correct theoretical model, even for high mass candidates. 
This is especially true for argon, which has the fewest number of events and smallest energy 
range. Xenon, however, has sufficient statistics to be able to extract the correct operator 
mediating the scattering on its own. The good discrimination of xenon is what makes the 
overlap between the xenon, germanium, and argon CLCs a sensitive test. 

We now investigate how CLC separation occurs. In Fig. ([T]) we compare input data corre- 
sponding to standard momentum independent scattering to scattering from q^'^ interactions. 
We specifically choose g^^ because these operators have the most exaggerated kinematic 
effects. These comparisons will lead to the largest separations in parameter space and show 
the widest disparity in the values of L^in/ d.o.f. The differential event rate of the mismatched 
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Figure 1. 95% CLCs for a 10, 50, and 250 GeV particle interacting through an rii = standard 
operator with equal exposures, such that there are 300 events on a germanium target and more 
(fewer) on a xenon (argon) target. Comparisons are made to Ui = standard, g^, and operators. 
The colors represent the value of Lmin/d.o.f. The standard interaction can be distinguished from 
the q^^ operators via both CLC overlap and the values of Lmin/d.o.f. As can be seen in Fig. (10) in 
the appendix, Lmin/d.o.f. is less powerful for distinguishing m = standard from Ui = g^^, though 
overlap remains robust. 
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Figure 2. 95% CLCs for a 10, 50, and 250 GeV particle interacting through an Ui = dipole operator 
with equal exposures. Comparisons are made to = standard, g^, and operators. The colors 
represent the value of -Lmin/ d.o.f. The dipole moment operator is dominated by a term proportional 
to g^, and is indistinguishable from this operator. The dipole interaction can be distinguished 
from the and standard operators via CLC overlap, which is largely powered by the statistics 
of the xenon target. As can be seen in Fig. (12) in the appendix, overlap remains powerful for 
distinguishing rij = dipole from other operators. 
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250 GeV, Argon spectra: n,— >std, v,— >o'* 250 GeV, Xenon spectra: «,-»std, Vi->q* 
7 40| | 
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Figure 3. Effect of momentum dependence and DM mass on fit spectra. The data are from an 
rii = standard interaction and are depicted with Poisson error bars. The best fit spectrum for 
Vi = standard is shown as a green sohd line. We also show two Vi = spectra. The spectrum 
in blue is at the best fit point and the q^ spectrum in red is taken at high mass to illustrate the 
improvement in goodness of fit from reducing uidm- We can also see the effect of xenon's higher 
event rate and lower energy threshold on its ability to determine the correct operator. 



operators is given by 

—— — [Er) ~ —= 5 — exp -7-^5 — 5- . (18) 

dER VvrmoMAi^t'o \ % J \ ^I^nK J 

The most noticeable effect of the additional momentum dependence, as shown in the second 
and third rows of panels of Fig. ([T]), is that on all targets the incorrect operators prefer 
different DM mass ranges than the value chosen by the correct operator: scattering 
chooses low mass and scattering prefers high mass. 

This behavior is driven by the shape of the event distributions. As seen in Fig. ([s]), 
the g"*"" operator with a lower mass candidate can mimic the scattering of a higher mass 
candidate with a standard spectrum. This occurs because as one decreases the event 
number increases at low energy and decreases at high energy. Likewise, we can consider a 
operator, shown in the bottom row of Fig. ([T]). In this case the event rate is suppressed 
at high recoil but increased at low recoil. In Fig. Q we see that larger DM masses are 
necessary to suppress the divergent low energy tail of the operator. 

Changing the momentum dependence by only two powers, e.g. comparing rii = std with 



15 



250 GeV, Argon spectra: n,— >std, Vi->q * 250 GeV, Xenon spectra: «;^std, Vi->q 
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Figure 4. Effect of momentum dependence and DM mass on fit spectra. The data are from an 
rii = standard interaction and are depicted with Poisson error bars. The best fit spectrum for t'j = 
standard is shown as a green solid hne. We also show two Vi = q~'^ spectra. The spectrum 
in blue is at the best fit point and the spectrum in red is taken at low mass to illustrate the 
improvement in goodness of fit from increasing rriDM- We can also see the effect of xenon's higher 
event rate and lower energy threshold on its ability to determine the correct operator. 



Ui = q , offers less contrast. In this case, fits for individual elements can in general be 
acceptable (with Lmin/d.o.f. < 1), and the overlap test seems to be an important discrim- 
inant. This can be seen for the = dipole case displayed in Fig. ([2]), where we compare 
the dipole input data to standard, g^, and interactions. We see that discrimination is 
marginal for an individual target (since Z/mm/d.o.f can be good), but the overlap appears to 
provide conclusive discrimination against the standard and g^-type interactions. By contrast, 
and dipole interactions look nearly identical. This is because the q^ term in Eq. ([T]) 
dominates the q^ term for the elements and DM masses examined here, as seen in Fig. ([s]). 
As expected, the standard spectra are discrepant with the dipole spectra due to a divergent 
low mass tail, the q'^ operator suffers from being overly suppressed at low energy, and the q^ 
operator provides very good fits. For the same reason, we see in the appendix that there is a 
degeneracy between the standard and anapole operators. Although the anapole operator has 
a non-standard velocity-dependence, its spectrum is very similar to momentum independent 
scattering. Since the velocity dependence does not affect the spectrum and the piece is 
subdominant, the anapole and standard operators may be said to have the same momentum 
dependence. Likewise, the dipole and g^ operators share their own momentum dependence. 

Based on the collected plots in the appendix, we can extend these conclusions. For 
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250 GeV, «,^dip, v/^std 250 GeV, «/^dip, Vi^q^ 250 GeV, «;^dip, v,^?* 
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Figure 5. Fits of Ui = std, g^, and operators to rii = dipole data. The data are shown as points 
with Poisson error bars. The dipole best fit spectrum is shown as a sohd hne and the other operator 
best fit spectra are shown as dotted hues. Xenon is red, germanium is blue, and argon is green. 



example, on argon, which has a high energy threshold and small event number, several 
incorrect operators can give a decent global fit to the data and a low value of Lmin/d.o.f. 
Even on germanium the wrong operator can mimic the true operator if the momentum 
dependence of the operator differs by a power of or less. This is not true for xenon, which 
has a good I/mm/d.o.f. only for the correct operator. Thus the Lmm/d.o.f. test works for 
xenon but not for argon. Additionally, it appears evident from the figures collected in the 
appendix that the overlap test, combining data from all targets, also works very well for 
high mass candidates. Nonetheless, the overlap test described here is in some sense another 
manifestation of the Lmin test. This comes about because of the relatively large number of 
events on the xenon target, which increase its value of Lmin (since, for large event numbers, 
Lmin scales with the number of events) and shrink its CLC relative to the other CLCs. The 
small xenon CLCs are in turn more precise in selecting a preferred region in parameter 
space, and thus are the single most important factor in the overlap test. In this way, the 
overlap and Lmm tests are linked. This can even be seen in the 10 GeV case, for example, 
by comparing the CLCs in the appendix from the = and standard data sets. In the 
former case the xenon cross section is enhanced relative to germanium by the ratio of the 
xenon and germanium atomic ratios Axe/^ce, which increases the xenon sensitivity, shrinks 
its CLCs, and allows modest CLC separation. In the rii = standard case there is no such 
enhancement and no CLC separation is evident. Now we turn to investigating the impact 
of statistical effects on our conclusions by considering the case where all targets observe the 
same number of events. 
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Figure 6. 95% CLCs for 10, 50, and 250 GeV particles interacting through an nj = standard 
operator with 300 events on each target. Comparisons are made to Vi = standard, q^, and 
operators. The colors represent the value of Lmin/ d.o.f. The standard interaction can be 
distinguished from the operators only via the values of Lmin/d.o.f., unlike in Fig. ([T]) where 
CLC overlap also provides a test. 
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Figure 7. 95% CLCs for a 10, 50, and 250 GeV particles interacting through an ni = dipole 
moment operator with 300 events on each target. Comparisons are made to = standard, g^, and 
operators. The colors represent the value of Lmin/d-o.f. Compared to the case shown in Fig. ^ 
where exposures were fixed for all targets so that xenon was able to power the statistics, the ability 
to discriminate the dipole interaction from standard and operators using overlap or the values 
of Lmin/d.o.f. is diminished. 
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C. Operator Discrimination with Equal Event Numbers on All Targets 

In Figs. ^ and ^ we show the CLC plots for 50 GeV and 250 GeV candidates with 300 
events on all targets to see if any discriminatory ability is lost. We omit the 10 GeV case 
because of the poor discrimination; we have checked that even by increasing the number 
of events in the 10 GeV candidate case by an order of magnitude, little improvement 
occurs. To achieve these event numbers we simply reduce the exposure for xenon; for argon 
we increase the exposure in addition to lowering the energy threshold from 20 keV to 5 
keV. This "equal event number" normalization is less physically motivated than the "equal 
exposure" normalization adopted above, but is necessary for understanding the robustness 
of our results. 

Compared to the figures from the equal exposure scenario it is apparent that for high 
mass candidates the simple overlap test loses much of its capability to distinguish operators. 
There is no overlap when Ui = standard and Ui = g^, but there is overlap for z/j = This 
problem is even more noticeable in the case with rij = dipole where with equal events there 
is mutual overlap in almost every instance. This loss of power is expected since the precision 
of the xenon CLCs has greatly diminished. 

However, the Lmin/d.o.f. test remains generally strong and even increases in relevance 
for some combinations of operators. Previously, only xenon and germanium were capable of 
differentiating operators on the basis of the Lmm/d.o.f. test. By increasing the number of 
events on an argon detector we have a third viable test of the Lmin/d.o.f. of each operator. 
This compensates for the weakening of the overlap test that is a consequence of decreasing 
the number of events on xenon. Yet, even though we are still able to extract the correct 
operator, it seems that our standards for discrimination have been changed by the statistics. 

We would like to be able to use the same test on any combination of operators, elements, 
and exposures. The fact that the overlap test is powered by the high Lmm values of xenon 
is an important clue. We will show that what seems to be a qualitative difference in how we 
draw conclusions in the equal exposure and equal event scenarios is just an illusion. 



D. Comparison of Models 



In this subsection, we present significance tests for different trial operators. We combine 
all three data sets and define i^mL be the minimum log-likelihood value from this global 
fit. Here we make the approximation that -Z^mL values follow a distribution. This is valid 
when there are a large number of events, as in the cases studied here. Thus, we may derive 
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a p- value for any given trial operator by defining [23] 

P'^"°-"-'' = 4-„ 2"^'rK/2) ''"- 

where is the number of degrees of freedom; this is the total number of nonzero bins for 
all three experiments minus the number of fit parameters. For a given trial operator, the 
p-value describes the probabihty of producing a fit that has a larger than L^ln, so that 
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Figure 8. The p-values of trial operators for n, = standard and operators for candidate masses 
of 50 and 250 GeV. We display the equal exposure and equal event bar charts side by side to 
underscore the robustness of the discrimination. For the visual purpose, the plot is normalized so 
that each bar starts at 5% significance level. 
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a high p-value for an operator indicates that the relevant model is able to fit the data well. 

We display results for the standard, anapole, dipole, g^, and data sets in Figs. ^ 
and ([9]). We omit the 10 GeV data from this comparison because the discrimination is 
obviously quite poor, and we are primarily interested in seeing if the higher mass candidates 
are sensitive to the statistical effect of changing the event numbers seen by the argon and 
xenon targets. One immediately notices that the scattering of 50 and 250 GeV candidates 
typically provides very good differentiation between correct and incorrect operators. In most 
of cases, we can reject the trial model at much better than 5% significance level if operators 
do not match. Also as expected, we see that the standard and anapole operators are hard 
to distinguish, as are the and dipole operators. This is a consequence of the fact that the 
velocity-dependent contribution to the composite operators' spectra has weak momentum 
dependence, and thus the anapole and standard operators have approximately momentum 
independent spectra just as the dipole and operators have the same momentum 
dependence. 

Modulo this degeneracy, we see that it is possible to extract the particle physics nature 
of the dark matter scattering events given that future ton-scale direct detection experiments 
observe (9(100) dark matter events. Moreover, the p-value is capable of extracting the 
momentum dependence of the operator equally well for the two event normalizations (equal 
exposures or equal event numbers on all targets) studied here. Finally, we also have checked 
that the results obtained from the p-value test derived from the global-fit L]^*^ agree with 
the conclusions we would draw from calculating the p-values for each element individually, 
but the result from the global-fit has a better discrimination ability. 



V. CONCLUSIONS 



We have discussed the capability of direct detection experiments to extract the particle 
physics underlying DM scattering events. We found that Poisson noise limits the ability 
of a single detector to determine the momentum dependence of the operator mediating the 
scattering. When the exposures are equal for different elements and the data are from 
observations of high mass candidates, only xenon, due to its high atomic number, is capable 
of determining the momentum dependence of the interaction on its own. Under these 
conditions, examining the preferred regions from data generated by more than one element 
allows one to extract the correct momentum dependence; this can be seen in Figs. ([T]) and 
When two operators give a similarly good fit to the data, it is because they have a 
very similar momentum dependence. We found that much of the power of this "overlap 
test" is derived from xenon's capability to discriminate operators. When all event numbers 
are held constant, so that the xenon exposure is decreased and the argon exposure is scaled 
up, we found that the overlap test is less powerful, but the combination of data from all 
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three targets still provides good operator discrimination. We found that the minimum log- 
likelihood value for the global fit summarizes all of this information succinctly and robustly. 
The minimum log-likelihood value for the global fit is sensitive both to overlap and to the 
individual goodness of fit information, and can be put in one-to-one correspondence with the 
p- value, allowing us to make quantitative statements about the ability of each operator to fit 
the data. For a light DM candidate on the other hand, the effects of the energy threshold, 
the energy resolution, and the scattering kinematics combine in such a way that the operator 
mediating the scattering cannot be extracted even for very high event numbers. 

We thank Tim Cohen for discussions. The work of SDM and KMZ was supported by NSF 
CAREER award PHY 1049896, and the work of HBY and KMZ was supported by NASA 
Astrophysics Theory Program grant NNX11AI17G. 
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Appendix A: Plots 

Here we compile plots for analyses of all five data sets (rij = standard, anapole, dipole, 
and with equal exposure. Due to space constraints, we show fits to six of the seven 
operators (z/j = standard, anapole, dipole, g^^ and q^^) on each page. More comprehensive 
plots, including plots from data sets with equal event numbers, are shown online |27] . 

The plots illustrate the qualitative and quantitative features described in Section IV 
above. Due primarily to the power of the discrimination of the xenon target, analysis that 
results in three overlapping CLCs will indicate an operator that fits the data well. For low 
mass DM all operators are effectively indistinguishable in the sense that all analyses result 
in overlapping CLCs. When the DM candidate is more massive, the distinguishing features 
of its interaction become more pronounced and harder to mimic, so a particular type of 
operator is selected. 

Within each plot we display CLCs for all three detector elements and all three DM 
candidate masses, with the exception of the argon contour with the 10 GeV candidate. None 
of the 10 GeV argon CLCs close because argon sees so few events. The argon "exclusion 
curve" fits comfortably around the xenon and germanium CLCs in all cases. To reduce 
clutter, these argon curves are omitted from the plots. 
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Figure 10. 95% CLCs for a 10, 50, and 250 GeV particle interacting through an ni = standard 
operator. Comparisons are made to Ui = standard, anapole, dipole, q^, g^, and operators. The 
colors represent the value of -Lmin/ d.o.f. As described above, cyan and lighter colors correspond to 
95% or worse disagreement with the data. 
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Figure 11. 95% CLCs for a 10, 50, and 250 GeV particle interacting through an Ui = anapole 
moment operator. Comparisons are made to = standard, anapole, dipole, q^, q^, and q~'^ 
operators. The colors represent the value of Lmin/ d.o.f. As described above, cyan and lighter 
colors correspond to 95% or worse disagreement with the data. 
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Figure 12. 95% CLCs for a 10, 50, and 250 GeV particle interacting through an = dipole moment 
operator. Comparisons are made to Vi = standard, anapole, dipole, g^, g^, and operators. The 
colors represent the value of -Lmin/ d.o.f. As described above, cyan and lighter colors correspond to 
95% or worse disagreement with the data. 
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Figure 13. 95% CLCs for a 10, 50, and 250 GeV particle interacting through an rii = operator. 
Comparisons are made to vi = standard, anapole, dipole, q^, g^, and operators. The colors 
represent the value of Lmin/ d.o.f. As described above, cyan and lighter colors correspond to 95% 
or worse disagreement with the data. 
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Figure 14. 95% CLCs for a 10, 50, and 250 GeV particle interacting through an m = operator. 
Comparisons are made to fj = standard, anapole, dipole, q^, and operators. The colors 
represent the value of Lmin/ d.o.f. As described above, cyan and lighter colors correspond to 95% 
or worse disagreement with the data. 



32 



